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Since the work of Anderson on localization, interference 
effects for the propagation of a wave in the presence of dis- 
order have been extensively studied, as exemplified in coher- 
ent backscattering (CBS) of light. In the multiple scattering 
of light by a disordered sample of thermal atoms, interfer- 
ence effects are usually washed out by the fast atomic mo- 
tion. This is no longer true for cold atoms where CBS has 
recently been observed. However, the internal structure of 
the atoms strongly influences the interference properties. In 
this paper, we consider light scattering by an atomic dipole 
transition with arbitrary degeneracy and study its impact on 
coherent backscattering. We show that the interference con- 
trast is strongly reduced. Assuming a uniform statistical dis- 
tribution over internal degrees of freedom, we compute ana- 
lytically the single and double scattering contributions to the 
intensity in the weak localization regime. The so-called lad- 
der and crossed diagrams are generalized to the case of atoms 
and permit to calculate enhancement factors and backscatter- 
ing intensity profiles for polarized light and any closed atomic 
dipole transition. 



I. INTRODUCTION 

Interference of waves is the general feature shared by 
different fields of physics such as Optics, Acoustics and 
Quantum Mechanics. For waves propagating in disor- 
dered media, it was believed that interference effects 
would be scrambled and that a reliable Boltzmann trans- 
port theory would emerge m . But Anderson [EJ predicted 
in the context of solid state physics that interference can 
inhibit the propagation of matter waves in disordered me- 
dia (Anderson localization) . Since then, many theoretical 
and experimental works have shown that elastic multi- 
ple scattering in the presence of disorder is full of rich 
phenomena &-01. The coherent backscattering effect, an 
interferential enhancement of the average reflected light 
intensity in the backscattering direction, was the first di- 
rect experimental evidence H-H] that interference of light 
waves persists in the presence of disorder and has been 
extensively studied for the past fifteen years. 

At the same time, considerable advances were achieved 
in creating and controlling dilute gases of cold atoms, 
leading to the experimental observation of Bose- Einstein 
condensation in 1995 and triggering active experimental 
and theoretical research [gj . It is not surprising that cold 
atomic gases have been suggested as promising media for 
strong (Anderson) localization of light [nO[ . Well-defined 



atomic transition lines allow strongly resonant light scat- 
tering with cross-sections of the order of the squared op- 
tical wavelength, much bigger than the actual size of the 
atom. In this respect, atoms are natural realizations of 
the mathematical concept of point dipole scatterers (also 
known as resonant Rayleigh scatterers), a paradigmatic 
model in the context of multiple scattering [|ll|--|l3| . How- 
ever, this simplified description has become questionable. 
The atomic dipole transition interacting with light in 
real experiments is usually more complicated: both the 
ground state (with angular momentum J) and the ex- 
cited state (with angular momentum Jo) present an im- 
portant degeneracy, necessary for cooling and trapping. 
This internal structure makes the atom behave very dif- 
ferently from a point dipole scatterer. Indeed, coherent 
backscattering of polarized light by a laser-cooled gas of 
Rubidium atoms has been observed recently |l3,n5| in 
the weak localization regime. There, surprisingly low 
enhancement factors for the backscattered intensity in- 
dicate that interference is less efficient for atoms than for 
classical point dipole scatterers. A careful study of the 
coherent propagation of light waves in atomic gases there- 
fore promises to be of great interest for both fields "mul- 
tiple scattering in disordered media" and "cold atoms" . 

In this paper, we show in detail how the internal atomic 
structure can account for the reduction of the enhanced 
backscattering of polarized light by atoms. In particu- 
lar, we generalize the theory of single and double scat- 
tering of polarized light by classical point scatterers to 
the case of atomic scatterers with an arbitrarily degen- 
erate dipole transition. Because of this degeneracy, the 
full atomic scattering tensor has to be considered. It will 
be shown that its non-scalar parts are responsible for a 
single scattering background in all polarization channels 
and a drastic reduction of the interference contrast. 

The paper is organized as follows. Sec. O introduces 
the basic notions of enhanced backscattering of light by a 
standard disordered medium. Sec. Ill presents an analy- 



sis of single and double scattering amplitudes of light by 
atoms and shows qualitatively how a quantum internal 
structure reduces the backscattering enhancement. In 
Sec. IV, the single and double scattering intensities are 
calculated analytically, preparing the way for the quan- 
titative analysis contained in Sec. |V|. 



II. ENHANCED BACKSCATTERING OF LIGHT 



B. Enhancement factor 



A. Two-wave interference 

A wave, characterized by its wavelength A or wavenum- 
ber k = 2tt/X in vacuum, incident upon a disordered 
medium is scattered into a multitude of partial waves. If 
the individual scattering events are elastic, these partial 
waves are all coherent and interfere. In the weak localiza- 
tion regime, individual scatterers with a scattering cross- 
section a are distributed with number density n so that 
the scattering mean free path i — 1/na is much larger 
than A. This condition, equivalently stated as fc£ ^ 1, 
says that the mean distance between scattering events is 
much larger than the wavelength, so that waves propa- 
gate almost freely inside the medium. In this regime, the 
wave amplitude A can be constructed by coherent super- 
position A = ^ ap of partial waves which are scattered 
along a quasi-classical path joining the positions of con- 
secutive scatterers [Q. The positions of all scatterers 
in turn determine the precise shape of the resulting in- 
terference pattern, as observed in the speckle figures of 
scattered laser light. 

This interference pattern is naively expected to be 
washed out when averaged over the realizations of the 
disorder (for example, by thermal motion of the scatter- 
ers). In fact, the average intensity / — (|^P) separates 
into independently squared amplitudes and the sum of 
interference terms, / = X]p (lopP) + Hp^p' (o-pap') (the 
brackets indicate an average over realizations of disorder, 
the bar denotes complex conjugation). If the scatter- 
ers are distributed randomly, different scattering paths 
(p' y^ p) involve uncorrelated phases. The interference 
terms may be expected to vanish, {apUp') = 0, yield- 
ing the uniform average intensity that is familiar to us 
from the view of most natural objects like clouds. In 
the context of light scattering, it was first realized by 
Watson [0, de Wolf Q and others, however, that each 
multiple scattering sequence visiting N scatterers in a 
given order (1, . . . , A^) has exactly one reversed counter- 
part {N, . . . , 1). The phase difference between the two 
corresponding partial waves (visiting the same scatter- 
ers, but travelling in opposite directions) is given by 
A(f) = (k -I- k') • (ri — r^r), where k and k' are the in- 
coming and outgoing wave vectors, and ri and tn are 
the positions of the first and last scatterer along the 
scattering path. The phase difference is exactly zero in 
the backscattering direction where k' = — k. Zero phase 
difference means constructive interference, independently 
of the actual path configuration. This constructive two- 
wave interference therefore survives the ensemble average 
and gives rise to coherent backscattering, the enhance- 
ment of the multiply scattered intensity in the backward 
direction by a factor of two. 



There is an exception to the systematic interference 
between direct and reverse amplitudes: scattering paths 
which are their own reversed do not give rise to any in- 
terference term and thus add a uniform background to 
the average scattered intensity. In the weak scattering 
regime ki ^ 1, this uniform background reduces to the 
single scattering contribution /g. In this regime, the av- 
erage intensity can be written as a sum of three terms 
I{e) = Is{e) + Il{9) + Ic{9) as a function of the an- 
gle 6 with respect to the backscattering direction. Here, 
the so-called ladder term /l {0) is the contribution of all 
squared multiple scattering amplitudes, neglecting inter- 
ference. The so-called crossed term /c(^) contains the 
interferences between direct and reverse amplitudes. Un- 
der well chosen experimental conditions, where all paths 
and their reverse counterparts have exactly the same am- 
plitude, the constructive two-wave interference leads to 
a maximal contrast /c(0) = /l(0). Away from the back- 
ward direction, Ic{d) is averaged to zero once the typi- 
cal phase difference of interfering amplitudes approaches 
A0 K, 1. Taking the double scattering contribution 
[N — 2), the distance ri2 = ||ri — r2|| will be given 
on average by the scattering mean free path L To first 
order in 0, the typical phase difference then is A(/) — k£9. 
Therefore, Ic{6) decreases to zero over an angular scale 
1/ki which is very small in the weak scattering regime 
fc£ 3> 1. Higher orders of scattering involve paths with 
endpoints further apart and thus contribute to Iq with 
a smaller angular width. For a semi-infinite and non- 
absorbing scattering medium, the sum of all contribu- 
tions has been shown to result in the so-called coherent 
backscattering cone, a sharp intensity peak exactly in the 
backscattering direction |19|-pi|. When higher orders of 
scattering become relevant, the width of the backscat- 
tering enhancement is determined not by the scattering 
mean free path but rather by the transport mean free 
path £ti = £/{! — (cos 6*)); here, (.) denotes an aver- 
age over the differential cross-section. If {cos 9) = 0, 
the two length scales are identical, Itr — ^- This is true 
for i sotrop ic point scatterers and unpolarized atoms (cf. 
Sec. IV D), so scattering and transport mean free path 
will be identified throughout the rest of this article. Is 
and /l exhibit a smooth angular dependence with re- 
spect to the normal of the surface of the medium (Lam- 
bert's law p2[). They can thus be taken constant, for 
not too oblique incidence, on the backscattering angular 
scale l/k£. 

The ratio of the average intensity at backscattering 
1(0) = Is+Il+Ic{0) to the average background intensity 
I{k£9 :^ 1) ~ Is + II ^s, the enhancement factor 

» = i + ^. (1) 

JS + ^L 

Its maximum value a = 2 is attained if and only if there is 
no single scattering background. Is = 0, and the contrast 
of the two- wave interference is perfect, /c(0) = II- 



C. Polarization and reciprocity 

Since light is a vector wave, polarization (which de- 
scribes the direction of the electric field vector) is an es- 
sential ingredient of any analysis of the enhancement fac- 
tor. The incident field polarization e and the scattered 
field polarization e' define two sets of orthogonal polar- 
ization channels. For linearly polarized incident light, the 
scattered light can be analyzed with parallel [lin || lin) or 
perpendicular {lin _L Uri) polarization. For circularly po- 
larized incident light, it is convenient to use the concept 
of helicity, i.e., the orientation of the circular polarization 
with respect to the direction of propagation. The scat- 
tered light can be analyzed with preserved helicity (ft, || K) 
or flipped helicity (hJ-h). At exact backscattering, these 
two cases respectively correspond to flipped (e' = e) and 
preserved polarization (e' = e). Note that the circularly 
polarized light scattered backwards by a mirror has the 
same polarization, thus flipped helicity. 

For classical scatterers, the following results have been 
established |2^,p3]: (i) Single scattering in the backscat- 
tering direction is absent in the lin _L lin and ft || ft, chan- 
nels for scatterers of spherical symmetry; (n) In the ab- 
sence of an external magnetic fleld, the reciprocity the- 
orem (see below) assures that /c(0) — /l in the parallel 
channels lin \\ lin and ft || ft. Satisfying simultaneously 
conditions (?) and (ii), an enhancement factor a = 2 has 
been predicted and observed for spherically symmetric 
scatterers in the ft || ft polarization channel pq]. 

As reciprocity is an important notion for CBS, let us 
precise this point. Reciprocity is a symmetry property 
stemming from the invariance of the fundamental mi- 
croscopic dynamics under time-reversal [ p4[ . Reciprocity 
assures that amplitudes relating to scattering processes 
where initial and final states are exchanged and time- 
reversed are equal. For the scattering of incident light 
with wavevector k and polarization e into light with 
wavevector k' and polarization e', it implies 



rdh.(k£ ^ k'e') = T„v(-k'£' ^ -k£) 



(2) 



Here, rdir(k£ -^ k'e') is the amplitude of a given scat- 
tering sequence, and Ti-ov(— k'£' -^ — k£) is the ampli- 
tude of the reciprocal process (the bar denotes complex 
conjugation). In general, these reciprocal amplitudes 
describe different scattering processes and thus cannot 
interfere. CBS interference arises between amplitudes 
'?dir,rov(k£ -^ k'e') associatcd to direct and reverse scat- 
tering paths with the same initial and final direction of 
propagation and the same polarization. The reciprocity 
relation (g) thus assures equality for the two CBS ampli- 
tudes if and only if two conditions are met: 



k = — k and 



£ = £. 



(3) 



^From these conditions, it follows that the CBS ampli- 
tudes of any given path and its reverse are equal at 
backscattering in the lin \\ lin and ft || ft channels, im- 
plying /c(0) = Ih- On the other hand, away from the 



backscattering direction or in the perpendicular chan- 
nels, the relation (||) is still valid, but says nothing about 
the pairs of amplitudes that interfere for CBS. These am- 
plitudes are therefore expected to be different, leading to 
a reduced contrast /c < II- 



III. AMPLITUDES FOR SCATTERING OF 
LIGHT BY ATOMS 

A. Description of the atomic medium and 
approximations 

We are interested in the situation where the scatter- 
ers are not macroscopic objects, but individual atoms. 
One may think of several specific characteristics of atomic 
light scatterers which affect coherent backscattering: 

• Atoms have extremely narrow resonances. Close 
to an atomic resonance, the light scattering cross- 
section is of the order of the square of the wave- 
length, much larger than the geometric cross- 
section of the atom. A dense cloud of atoms there- 
fore is ideal for strong elastic multiple scattering. 

• Because of the high polarizability of atoms near 
an atomic resonance, it is rather easy to induce 
non-linear effects (e.g. saturation) with only few 
milliwatts of laser power. Despite some studies of 
multiple scattering in non- linear media p6| , it is 
basically unknown how this affects CBS by atoms. 

• When an atom scatters a photon, its velocity 
changes by an amount of the order of mm/s. This 
recoil effect becomes important for cold atoms with 
typical velocities of a few cm/s. 

• The atomic resonances being very narrow, atoms 
may be driven in or out of resonance because of 
the Doppler effect. Adding the contributions of the 
various velocity classes to the CBS signal is far from 
obvious. 

• Atoms also have a quantum internal structure. For 
a given transition line, the total angular momen- 
tum J of the atomic groundstate in general is not 
zero. In the absence of any external magnetic field, 
the groundstate then is (2J -I- l)-fold degenerate. 
As a first consequence, there is the possibility of 
elastic light scattering processes which change the 
internal atomic substate (degenerate Raman tran- 
sitions) . Subsequent light scattering then gives rise 
to optical pumping. 

• When the atoms are very cold, their de Broglie 
wavelength becomes comparable to the optical 
wavelength. In this regime, the external atomic 
motion must be treated quantum mechanically. For 
high enough density, Bose-Einstein condensation 
sets in. 



Addressing all these problems is beyond the scope of 
this paper. We will focus our present investigation on 
the crucial role of the atomic internal structure, making 
use of several simplifying approximations. 

Firstly, we assume the weak scattering relation k£ ^ 1 
to hold. This will be the case for sufficiently low density 
n of the atomic medium. Indeed, as the resonant atomic 
cross-section a — 1/nl is of the order of A^, weak scat- 
tering is implied by the low density condition nX'^ <^ 1. 
In this regime, the independent scattering approximation 
(ISA) is justified pl. Equivalently, recurrent scattering 
(i.e., sequences visiting a given scatterer more than once) 
can be neglected. The single scattering transition matrix 
then suffices to compute the single scattering intensity 
which, in turn, serves as a building block for higher order 
scattering. In this regime, the average index o f refra ction 
of the medium is very close to unity (cf. Sec. [V A ). 

Secondly, we use quantum mechanical perturbation 
theory to describe the scattering of light by an atom. 
This will be valid as long as the laser intensity is suffi- 
ciently low H^. We will restrict our calculation to the 
case of one photon scattering, determining the transient 
reponse of the system rather than its stationary state. 
This method thus ignores saturation effects and optical 
pumping. In principle both could be described by carry- 
ing the perturbation to higher numbers of scattered pho- 
tons, but practically one has to calculate the stationary 
density matrix by solving the multi-level optical Bloch 
equations. In the experimental application sofar p5| , the 
laser intensity was kept well below the saturation inten- 
sity. Furthermore, optical pumping in the bulk of an 
optically thick atomic cloud is expected to be severely 
limited by multiple scattering. 

Thirdly, we treat the external motion of the atoms 
classically. In other words, we require the atoms to be 
sufficiently hot so that the coherence length of the exter- 
nal wavefunction is shorter than the optical wavelength. 
This is the case for cold atoms created in a standard 
magneto-optical trap. The present treatment does not 
apply, however, to ultra-cold atoms as, e.g., in a Bose- 
Einstein condensate. 

The question of the recoil effect can then be addressed 
rather easily. Indeed, light imparts various momentum 
kicks to the atoms defining a scattering path, but these 
momentum transfers are identical for the direct and re- 
verse paths at backscattering. Consequently, the recoil 
effect does not affect the interference between the ampli- 
tudes along the two paths. 

Let us now discuss the role of the atomic motion. Close 
to an atomic resonance of width F, the average time an 
atom takes to scatter a photon is F^^. If, in the mean- 
time, atoms move by more than an optical wavelength, 
then the interference term between direct and reverse 
scattering sequences will be spoiled [g8|. To avoid this, 
we require the spread v of the atomic velocity distribu- 
tion to satisfy 



If this condition is met, atoms can be thought as being 
fixed in space during the multiple scattering process [[l5| . 
On a much longer time scale, the motion of the atoms 
simply acts as a configuration average. Typically, eq. (jj) 
is satisfied for atoms slower than few m/s, which is true 
for atoms originating from a magneto-optical trap. Note 
that eq. (m can be alternatively viewed as a resonance 
condition: under scattering, the Doppler shift will not 
bring atoms out of resonance. 

Under these conditions, the most important effect will 
come from the internal structure of the atoms, i.e., the 
degeneracy of the light scattering transition. We as- 
sume that the incident light is nearly resonant with an 
atomic transition of (bare) angular frequency ujq between 
a ground state with total angular momentum J and 
an excited state with total angular momentum J^ (see 
fig. |lj). Since no external magnetic field is supposed to 
be present, the atomic ground state and the excited state 
are respectively (2 J -I- l)-fold and (2 Jc -I- l)-fold degen- 
erate. The corresponding substates with magnetic quan- 
tum numbers m and rUe with respect to some arbitrary 
quantization axis are denoted by \Jm) for the ground 
state {—J<m<J) and by | Jc'tIc) for the excited state 
(-Je <rne< Jc). 

The restriction to a single J — » Jc transition could 
be relaxed at the price of more complicated calculations 
since, in essence, the various transitions contribute inde- 
pendently to the atomic scattering tensor, the essential 
ingredient of our analysis as shown below. We will also 
assume that the J —> Je transition is closed, so that light 
scattering is purely elastic. Again, different final states 
with different energies could be included along the same 
lines of reasoning. 



Jeirie) 



UJo 




kv < F. 
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FIG. 1. Energy representation of a degenerate atomic 
dipole transition, here for J = 1, Jo = 2. Arrows mark atomic 
transitions from the initial substate, here \J,m = J), under 
scattering of a photon. In the absence of an external magnetic 
field, all transitions are elastic. Solid arrows: Rayleigh tran- 
sitions, conserving the magnetic quantum number (m' = m). 
Dotted arrows: degenerate Raman transitions, changing the 
magnetic quantum number (here m' = m — 1 and m' = m~2). 

In the following, we recall the amplitude for the scat- 
tering of one photon by one atom (single scattering) and 
determine the amplitude for the scattering by two atoms 
(double scattering). We discuss how the degeneracy of 
the atomic dipole transition affects the CBS enhance- 
ment factor (^. We will use a full quantum mechani- 
cal treatment of both atoms and electromagnetic field. 
While the internal atomic degrees of freedom can only 



be described quantum mechanically, the electromagnetic 
field will be described by quantum Fock states for rea- 
sons of symmetry. An equivalent treatment can be set 
up for low-intensity coherent states which are known to 
correspond closely to a classical light field p3 . 

Throughout the paper, transitions between identical 
atomic substates (m' — m) are called Rayleigh transitions 
and transitions between different substates (m' ^ m) 
are called degenerate Raman transitions. Let us stress 
that, since one-photon scattering on a degenerate atomic 
dipole transition is necessarily elastic, inelastic processes 
(also known as Raman scattering) are completely ab- 
sent of our analysis. In the following, we use natural 
units where h = c = 1 so that [length] = [time] = 
[frequency] ^ ^ = [energy] ^^. 



B. Single scattering amplitude 

In the single scattering situation, an atom at fixed po- 
sition r is exposed to a plane light wave with wave vector 
k, angular frequency uj — k and transverse polarization 
£. We describe the uncoupled system by the sum of the 
atomic internal hamiltonian and the free field hamilto- 
nian. 



i/o = ^0 7 I Jc771o)(Jomc| + 



k,E_Lk 



Jke^ke- 



(5) 



Here, Oke and a\^^ are the usual annihilation and creation 
operator of a transverse field mode with wave vector k 
and polarization vector e. The corresponding one-photon 
Fock state will be denoted |ke) where the transversality 
(k • s) = is understood. The interaction between atom 
and light is given in the dipole form by ^ = — D • E(r). 
The atomic dipole operator D connects the subspaces 
Tij and Tij^ (since we consider a closed transition, no 
other subspaces are involved) with reduced matrix ele- 
ment (Jo||D||J> = D^/2Jc + 1 [|o). The electric field 
operator at the atomic position is given by 



E(r) = i ^ f^Skeake exp[i(k • r) 

k,e±k 



The field strength E^j = (cj/2eoi'^)^/^ is defined in terms 
of a quantization volume L^ that eventually disappears 
in results of physical significance. 

The probability amplitude for a transition from an ini- 
tial state \i) = [ke; Jm) to a final state [/) = [k'e'; Jm') 
is the element Sfi of the scattering matrix. The transi- 
tion amplitude for i ^ f is written 



Sfi = —2iTTS{Lu — uj')Tfi{uj + iO), 



(7) 



in terms of the transition operator T. Here, because 
the atomic ground state is degenerate, energy conserva- 
tion, assured by the delta distribution, implies elastic 
light scattering (to' — to). The matrix element Tfi{z) 



is calculated perturbatively using the Born expansion 
T{z) = V + VGo{z)V + . . . in powers of the interaction V 
and the resolvent Go{z) = {z — Hq)^^ of the unperturbed 
system. The excited atomic state can be eliminated by 
partial summation of the Born series, dressing the tran- 
sition frequency and introducing a finite lifetime [p7[ : 
(5 = w — Wo is the detuning from the (dressed) transition 
frequency uJq « wq and F = D^luq/Stteo is the natural 
width of the atomic excited state. 

Let us define the reduced dipole operator d = D/D 
and the Rabi frequency wr = DS^^. The transition ma- 
trix element Tfi = (k'e', Jm'\T{uj + iQ)\\i.e, Jm) near res- 
onance then is 



Tf^ = 



6 + iT/2 



{Jm'\{e' ■ d)(£ • d)| Jm) e'^^-^">-\ (8) 



represented by its Feynman diagram in fig. y. The condi- 
tion "near resonance" means 6 ^ ujq (but not necessar- 
ily 5 < r). Therefore, anti-resonant scattering, i.e., first 
emission then absorption, can be neglected. 
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FIG. 2. Feynman diagram for the transition matrix ele- 
ment Tfi (eq. (0)) near resonance. Wavy lines denote pho- 
tons, straight lines atomic internal states. The thick line 
stands for the dressed propagator {S + iV/2)~^ of the ex- 
cited atomic state. The absorption vertex yields a factor 
cjR exp(ik • r)(Jeme|e ■ dj Jm), the emission vertex yields a 
factor LjR exp(— ik' • r){Jm'\e' ■ dj Jeme), and all intermediate 
variables, here ttIo, have to be summed over. 

In eq. (0), all information about the atomic internal 
degrees of freedom and polarization has been factorized 
into the matrix element 



(6) (Jm'|e'-t-£|Jm) = 



iF/2 



(Jm'|(r'-d)(e-d)|Jm). 
(9) 



This defines the scattering operator t which acts on the 
product space TLj ® C'^ of atomic internal states and po- 
larizations. It is the scattering operator t that character- 
izes the scattering object and contains all relevant infor- 
mation about the scattering process Q. We can sepa- 
rate its frequency dependence from its tensor structure: 
t{uj) = t{uj)t where 



t{uj) 



Wt 



5 + iT/2 



(10) 



is given as the ratio of the squared Rabi frequency (or 
squared coupling strength) and the resonant denominator 



known for point dipole scatterers ||ll|| . The novelty of the 
present approach hes in the pecuHar tensor part t. For 
a given transition tti — > m', the dimensionless matrix 
element 



iij{m,m') — {Jm'\didj\Jm) 



(11) 

defines the scattering tensor which connects the incoming 
to the outgoing polarization. This 3x3 i-matrix can be 
decomposed into its scalar, antisymmetric and symmet- 
ric traceless components, transforming irreducibly under 
rotations pOJ . 

A classical point dipole scatterer is characterized by 
a t-matrix proportional to unity |ll| . This behaviour 
is reproduced by the elementary dipole transition J — 0, 
Jo = 1. Indeed, the only matrix element of the scattering 
operator yields the scalar part (00|ty|00) = Sij. Non- 
spherical classical scatterers also display an additional 
traceless symmetric part in their scattering i-matrix. In 
the case of atoms, therefore, it is the antisymmetric part 
that is characteristic for the quantum internal structure. 
The antisymmetric part simply implies that an atom 
scatters light with polarization-dependent strength. To 
be more specific, consider scattering of circularly polar- 
ized light in a Rayleigh transition (m' = m; quantiza- 
tion axis is the direction of propagation) . The scattering 
strengths for the two possible helicities are different be- 
cause the Clebsch-Gordan coefficients associated to the 
transitions \J,m) <-^ \Jc,m± 1) are unequal. 

This situation is somewhat similar to the usual Fara- 
day effect where circular polarizations with opposite he- 
licities are scattered differently in the presence of an ap- 
plied magnetic field [^. There are, however, significant 
differences: in the Faraday effect, the antisymmetric part 
of the atomic polarizability depends both on the mag- 
netic field direction and on the direction of light propa- 
gation; for the atomic scattering operator, antisymmetry 
is a fully intrinsic property. When averaged over the in- 
ternal state, the antisymmetric part of the atomic scat- 
tering operator vanishes, leading to a symmetric polariz- 
ability a nd to no dichroism inside the effective medium 
(cf. Sec. IV A ). Thus, the degenerate atomic situation 
ressembles a (zero magnetic field) Faraday effect depend- 
ing on the internal state of the atom. 

The degeneracy of the atomic ground state also im- 
plies that, by scattering a photon, the internal state may 
change (cf . Fig. 1^) . The possibility of changing the inter- 
nal state leaves more choice for the photon polarization. 
Which polarization is possible for which transition follows 
from the conservation of angular momentum. In exactly 
the backscattering direction (k' = — k) and chosing the 
quantization axis along the direction of propagation, the 
following relations hold: linearly polarized light is scat- 
tered into the lin \\ lin channel by a Rayleigh transition 
(m' = m), and into the lin _L lin channel by a degenerate 
Raman transition {\m' — m\ = 1); circularly polarized 
light is backscattered into the hJ-h channel by Rayleigh 
transition (m' = m), and into the h\\h channel by a 
degenerate Raman transition {\m' — m\ = 2). 



Therefore, the single scattering amplitude shows that 
changes in the atomic internal state permit changes in 
the light polarization. Since in general the atomic inter- 
nal state is not under control, the single backscattering 
contribution cannot be removed by polarization analysis 
(with the only exception J = i in the h\\h channel) and 
degrades the observable enhancement factor (fil). 



C. Double scattering amplitudes 

1. Direct and reverse transition amplitudes 

In the double scattering situation, a plane wave im- 
pinges upon two atoms a = 1,2 at fixed positions r^. 
The interaction between atoms and field in the full haniil- 
tonian is now V — — Di • E(ri) — D2 • E(r2). This in- 
teraction defines a transition operator T along the lines 
of Sec. IIIB. Resonant dipole interaction between the 



atoms arises from the exchange of photons. Among the 
numerous different diagrams that describe the transition 
\i) — |k£, Jmi, JTO2) — > I/) = \]i'e',Jm[,Jm'2), the two 
dominant diagrams involving both atoms are concate- 
nations of two single scattering diagrams. The first di- 
agram, shown in fig. ^ describes the direct scattering 
path: absorption of the incident photon by atom 1 and 
emission of the final photon by atom 2. The diagram for 
the reversed path is obtained by exchanging the role of 
the two atoms. 



J nil 




Jm2 



FIG. 3. Feynman diagram of the direct transition matrix 
element Ti-" (eq. (|l2[)): resonant scattering first by atom 1, 
then by atom 2. The Feynman rules are defined in fig. ti. 

The Feynman rules introduced in fig. permit to ex- 
plicit the scattering amplitudes. As usual in diagram- 
matic expansions, a sum over the virtual intermediate 
states has to be carried out, here the excited atomic 
states and the intermediate photon. The matrix element 
for the direct scattering path in the far field approxima- 
tion kri2 ^ 1 takes the following form: 






3Ft2(a;) exp(ifcri2) 



4c^i, 



kri2 



idir e 



i(k-ri— k'-r2) 



(12) 



Again, all information about polarization and internal 
structure is factorizcd into the dimensionless matrix ele- 
ment 



idir = £' • {2(^2,^2) • A • ti(mi,m'i) • e. 



(13) 



Here, the dimensionless one-atom t-matrices, defined in 
eq. (|ll|), are connected by the projector Ay = 5ij — fiifij 
onto the plane transverse to the unit vector fi — ri2/ri2 
joining the two atoms. 

The matrix element T;'"^ for the reversed path is ob- 
tained from eq. (jl^) by exchanging the roles of atoms 1 
and 2. The internal matrix element ([l3[) becomes 



and a third one pertaining to the atomic internal vari- 
ables, 



ircv = e' • ti{mi,m[) ■ A • £2(7712, "^-2) ' £• 



(14) 



The two amplitudes Tf^"' and T^^^' describe indis- 
tinguishable processes and interfere. A maximal con- 
trast in the backscattering direction is obtained if and 
only if the amplitudes have equal magnitude idir — ^rcv. 
But due to the non-scalar part of the atomic i-matrix, 
we expect that in general the matrices do not commute, 
£2 • A • ti 7^ ti • A • t2, so that (^ and ([l|) are not 
equal. An exception to this rule is of course the case 
nil = 7ft2, iTT-'i = 7772 whcrc the exchange symmetry as- 
sures their equality. Furthermore, we can see that it is 
precisely the antisymmetric part of the i-matrix that is 
responsible for the inquality of amplitudes in the parallel 
polarization channels. Indeed, if the one-atom i-matrix 
were symmetric, then e' • ti(7r7i, m'l) ■ A • £2(7772, 7772) ■ £ — 
£ ■ t2 (7772, 7712) • A • ti (r77i , m'l) ■£'. In the parallel channels 
e' = £, from this would immediately follow the equality 



of the direct and reverse matrix elements ([TSf) and (14) 



But because of the antisymmetric part of the atomic scat- 
tering tensor, in general 



e' • t\{m\,m!i) ■ A • £2(7772,7712) • £ 7^ 



£ • £2(7772,7772) • A -£1(7771,7)7'^) • r'. 



(15) 



so that the two interfering amplitudes are different in 
magnitude. An explicit example for unequal interfering 
amplitudes in the h\\h channel - one is zero while the 
other is not - has been given in [p3|. 



2. Reciprocity revisited 

A question may arise at this point: Does the imbalance 
of amplitudes idir 7^ ^rov contradict the theorem of reci- 
procity? The answer is no: the complete system "field 
and atoms" obeys reciprocity, but this does not imply 
^dir — irev The classical reciprocity relation (0) has to 
be generalized to take into account the set {777} of the 
internal variables of all atoms |34| : 







rdU.(k£,{777} ^ k'£',{777'}) = (-1)^^^',' 

xT„v(-k'£', -{m'} -^ -k£, -{777}). 



(16) 



This relation shows that in order to obtain the recipro- 
cal sequence of a given sequence, the signs of all internal 
quantum numbers have to be flipped. The reciprocity 
relation (16) assures the equality of two interfering CBS 
amplitudes only if three conditions are met: the two clas- 
sical conditions (pi) on light direction and polarization. 



{m'} = {-777}. 



(17) 



Whereas the direction of observation and polarization 
can be controlled experimentally, this is impossible for 
the internal atomic states in an optically dense medium. 
Just as in the case of scattering away from the back- 
ward direction or into perpendicular polarization chan- 
nels, reciprocity including the internal states never ceases 
to be valid, but simply becomes inapplicable to predict 
the equality of interfering amplitudes. It follows that al- 
though there might be some amplitudes satisfying condi- 
tion (|l^), the majority will not, and perfect interference 
contrast is lost. Of course, in the case of the elementary 
dipole transition J — 0, Jc — 1, the condition (|l7| ) is 
trivially fulfilled since all atoms verify m' = m — and 
we recover the classical case. 

Classical reciprocity for light scattering has been de- 
rived using Maxwell's equations for a linear scattering 
medium provided that its constitutive tensors (dielec- 
tric constant, permeablility and conductivity) be sym- 
metric p5| . In the present case, when one does not con- 
sider the internal atomic states as intrinsic variables of 
the system but as given parameters for each path, the 
i-matrix £(777, 777') then has an antisymmetric part. In 
this respect, atoms with degenerate transitions consti- 
tute a scattering medium that does not obey classical 
reciprocity, and a reduced interference is no surprise. In- 
deed, the same is observed in scattering media with the 
Faraday effect |36y37|] where the external magnetic field 
is said to break time-reversal invariance. 

Finally, let us note that the ensemble average over the 
internal variables {777} cannot restore the equality of the 
direct and interference contributions to the diffuse inten- 
sity. Indeed, /c(0) is equal to Jl if and only if, for each 
pair of scattering paths, the direct and reverse amplitudes 
are equal. In the sum of all contributions, the equal am- 
plitudes cannot win back what the non-equal amplitudes 
have lost: the result is a reduced overall enhancement 
factor. 



3. The role of degenerate Raman transitions 

In the case of the elementary dipole transition J = 0, 
Jc = 1, the atomic scattering tensor only has a scalar 
part (00|iy|00) — 5ij. T he anal ysis of the double scat- 
tering amplitudes in Sec. [II C 1 shows that the internal 
amplitudes then are equal and full interference contrast 
is guaranteed. As for J = no degenerate Raman tran- 
sitions (777' 7^ m) can occur, the following question is 
inevitable: can the decrease of interference contrast be 
attributed to the Raman transitions alone? 

Indeed, one might be tempted to suggest incoherence 
of the Raman scattered light (i.e., the loss of phase co- 
herence in spontaneous emission) as the origin for this 



loss of contrast. In the present description, however, this 
is not a pertinent explanation. It is true that the Raman 
scattered light does not interfere with the reference light 
from the source: the respective final atomic states are 
orthogonal and the two amplitudes do not describe indis- 
tinguishable processes (this is a typical "which-path" ar- 
gument P8| , p9[ ). But a photon scattered elastically along 
the direct path interferes very well with the same photon 
scattered along the reverse path — as long as the inter- 
nal states of all atoms in both processes are identical, 
no matter whether they describe degenerate Raman or 
Rayleigh transitions. 

A closer analysis of the situation in the channels of cir- 
cular polarization permits the following remarks. In the 
channel /i _L /i of flipped helicity, a selection rule special to 
the double scattering configuration admits only Rayleigh 
transitions to the crossed intensity. The ladder intensity 
contains a contribution from Rayleigh transitions (equal 
to the crossed intensity) and an additional contribution 
from degenerate Raman transitions. In this sense, the de- 
generate Raman transitions are responsible for a reduced 
double scattering interference in the h J-h channel. This 
is consistent with the observation that degenerate Ra- 
man transitions make the atoms behave as non-spherical 
scatterers for which reduced interference in the perpen- 
dicular channels is expected. But for higher scattering 
orders, Raman transitions contribute also to the crossed 
intensity, and it is no longer evident to compare the rel- 
ative weights of Rayleigh and Raman contributions. 

The explicit example of a double Rayleigh transition in 
the h\\h channel with zero interference given in ref. [ p3[ , 
shows that Rayleigh transitions also are responsible for 
a loss of contrast. On the other hand, the fact that 
Raman transitions give to atoms some characteristics 
of non-spherical scatterers does not by itself imply a 
loss of contrast: the reciprocity theorem is indepen- 
dent of the actual shape of the scatterers and applies 
to spherical as well as to non-spherical classical scatter- 
ers. For example, a double Raman transition such that 
{mi = m2 = —m'l = —m'2 ^ 0) satisfies the reciprocity 
condition (^^ and has perfect contrast (as is also evi- 
dent from the exchange symmetry). In the sum of all 
scattering amplitudes, Raman scattering amplitudes can 
even be dominant in the backscattering interference sig- 
nal. An explicit example of such a situation is given in 

fig- I 

In fact, independently of the scattering order, it is pre- 
cisely the antisymmetric part of the atomic scattering 
tensor that is responsible for the loss of contrast in the 
parallel polarization channels. This antisymmetric part 
appears for both degenerate Raman and Rayleigh tran- 
sitions (cf. the unequal scattering of circularly polarized 
light with different helicities as mentioned in Sec. IIIB| ). 
Therefore, the degenerate Raman transitions must not 
be held responsible alone for the reduction of interfer- 
ence contrast. 



IV. ANALYTICAL FORMULATION OF THE 
INTERNAL ENSEMBLE AVERAGE 

We wish to describe the light propagation inside a 
macroscopic disordered medium on average. Starting 
from an entirely symmetric microscopic description of 
matter and light, the ensemble average is a trace over 
the matter degrees of freedom. This trace contains an 
average over atomic positions as well as an average over 
the internal degrees of freedom. This is analogous to the 
case of classical non-spherical scatterers where averages 
over position and orientation have to be performed. We 
suppose in the following that the atomic sample is pre- 
pared without correlations between positions and inter- 
nal substates, and that different atoms are uncorrelated. 
This is a reasonable assumption for a cloud of cold atoms 
created from a standard magneto-optical trap. The two 
averaging procedures then become independent. As far 
as positions are concerned, we will use the averaging tech- 
niques developed for classical point scatterers B . In the 
independent scattering approximation, the average over 
the internal quantum numbers {m} can be expressed as 
traces over a one-atom density matrix p and one-atom 
operators. 



A. Average amplitudes: effective medium 

Tracing over the matter degrees of freedom defines an 
effective medium for the average propagation of light am- 
plitudes and intensities. In this paragraph, we will deal 
with the rather simple issue of the average amplitude. As 
will be seen, the internal structure of the atomic scatter- 
ers provides no major surprise, and we are able to recover 
the well-known properties of a dilute atomic gas [ p9[ . 
The impact of the effective medium on the amplitude is 
described by the self-energy 'Sjuj) that renormalizes the 
vacuum light frequency lu |11| . In the independent scat- 
tering approximation, the self-energy is proportional to 
the average scattering operator. 



E(c.)-A^Trpt(c.)=iV(t(c.));„ 



(18) 



Because of the vector character of the light wave, the 
self-energy is formally a second rank tensor. Assuming 
a scalar density matrix, i.e., a uniform distribution over 
internal states, the internal average simply projects onto 
the scalar part: S(a;) = S(a;)l. Calculating the aver- 
age is elementary using the closure relation of Clebsch- 
Gordan coefficients, and we find 



E(w) = nMj 



Stt r/2 



(19) 



uj'^6 + iT/2 
Here, we define for convenience the ratio of multiplicities 



Ah 



2.L 



1 



3(2J-f 1) 



(20) 



with Mq — 1. The atomic polarizabihty close to reso- 
nance is given by 



a{uj) 



2L^ 



(t(^))i 



(21) 



Exphciting the internal average as a weighted sum over 
substates, 



(t('^)>mt =I]P" 



,(Jm|t(a;)| Jto), 



(22) 



it is evident that solely the Rayleigh transitions (m' = 
m) enter into the definition of the self-energy and of the 
polarizabihty. In the case of a uniform distribution with 
weights Pm = (2 J -I- 1)"^, this average selects the scalar 
part of the scattering operator. A thorough discussion of 
the polarizabihty, the scattering operator and its analysis 
by decomposition in irreducible components can be found 
in the textbook by Bcrcstctskii ct al. pl| . 

The susceptibility of the dilute atomic medium is 
X — na, and the condition of low density now reads 
na <C 1. The effective refractive index then is given 
by nr = 1 + na/2. Its real part is very close to unity, and 
we need not distinguish between the optical wavelength 
in the medium and in the vacuum. Its imaginary part 
describes attenuation of the average amplitude, and here 
the effect of the dilute medium is essential. Since we do 
not describe any absorption, all attenuation is necessarily 
due to scattering from the initial mode into other modes. 
This argument is the essence of the optical theorem 



tTtot = -2L^Im (e • t(w) • e)-^^ = k Ima(w) 
We thus find the total scattering cross-section 



CTtot = Mj 



Gtt 



1 



fc2 1 + A5^/V^ 



(23) 



(24) 



This well known expression features the resonant dipole 
cross-section Qir/k'^ — 3X^/2tt and the Lorentzian line 
shape for detuning 6 around the resonance with width F. 
The mean free path of the light inside the average 
medium is ^ = — 2(ImE(w))~^. By virtue of the opti- 
cal theorem, it depends on the total cross-section and on 
the number density of scatterers through: 



1 



natot 



(25) 



and is independent of both the polarization and the di- 
rection of propagation. This reflects statistical invariance 
of the atomic medium under rotation. 

In summary, in the weak density and weak scattering 
regime, the internal structure has very small influence on 
the properties of an average light amplitude. For a uni- 
form statistical distibution over internal states, all aver- 
age quantities are isotropic and are only modified by a 
factor Mj = (2 Je -I- 1)/3(2 J-|- 1) with respect to the clas- 
sical dipole point scatterer where Mq — 1. This is not 



surprising since the internal average over a scalar density 
matrix simply selects the scalar part of the atomic scat- 
tering operator. The antisymmetric part that appeared 
as the genuine quantum feature in Sec. Ill therefore is 
not present here. 



B. Average intensities 

Coherent backscattering is an interference effect for the 
average intensity, which of course must be distinguished 
from the square of the average amplitude. Consequently, 
we stress that it is not sufficient to calculate quantities 
pertaining to the average amplitude (such as the polar- 
izabihty or the scattering mean free path) in order to 
decide whether the internal structure affects CBS or lo- 
calization. In the following paragraphs, we show how 
proper use of tensor algebra makes it possible to ana- 
lytically perform the averaging over the atomic internal 
degrees of freedom. In the specific case of a semi-infinite 
medium, exact averaging over the external position of the 
atoms is also possible for the single and double scattering 
contributions. 

The average scattered intensity / far away from the 
medium can be calculated in terms of the dimensionless 
bistatic coefficient [Eo| 



7 



An 
1 



da 

In 



(ke ^ k'e') 



(26) 



Here, the light incidence is supposed to be perpendicular 
to the surface A of the medium which will be taken to 
become the half space z > as A ^ oo. The average 
differential cross section is determined by Fermi's golden 
rule which reads 



da 



k'e') 



4tt' 



(|r(k'£',k£)p). (27) 



The square of the transition operator means explicitly 
|T(k'e',ke)|2 = {ke\Tiuj + iOy\k'e'){k'e'\T{uj + iO)\ke) 
and acts on the atomic states only. Note that the fac- 
tor L^ cancels with the inverse factor coming from the 
squared transition operator, so that the quantization vol- 
ume finally disappears. 

For single scattering, eqs. (0) and (m show that the 
internal average has to be taken over the square of the 
(dimensionless) scattering operator: 

(|e' • t . ep)^^^ = Tr[p(£ . d)(£' • d)(e' • d)(£ • d)]. (28) 

It is crucial that the average be taken over the square of 
the scattering tensor. This is not equivalent to taking the 
square of the average which is essentially the polarizabil- 
ity (pl|). Again, the trace over a scalar density matrix will 
select the scalar part of the averaged operator. But, as 
becomes evident in App. ^, now the antisymmetric and 
symmetric traceless parts can combine with their coun- 
terpart in the direct product and contribute a non-trivial 
scalar component. 



In the double scattering situation, the two atoms are 
coupled by the intermediate photon. Let 



i2i = e^ • t2 • A • ti • £ 



(29) 



be a short-hand notation for the contracted double scat- 
tering operator for the direct path, and ii2 for the reverse 
path. The ladder contribution to the double scattering 
intensity, just like in the classical case, is given by the 
average sum of the squares of the two amplitudes, 



Tr[pi2(|i: 



211 



\ii2n. 



(30) 



Here, pi2 is the two-scatterer density matrix. The crossed 
contribution is obtained, again in perfect analogy to the 
classical case, by the interference between the direct and 
reverse amplitude 

Tr[pi2(ti2i2i exp[z(k + k') • ^2] + (1 ^ 2))]. (31) 

Since the atoms are uncorrelated, the density matrix fac- 
torizes, pi2 = pi ® P2- Furthermore, pa = p since the 
atoms are identically distributed. The two-atom scatter- 
ing operator (^9|) is not factorized in terms of elementary 
scalar products. But by expliciting the transverse pro- 
jector Ay — Sij ~ fiifij^ it becomes 

t2i - (£' • d2)[(d2 • di) - (n • d2)(n • di)](£ • di). (32) 



All averages ( pSD , (|30[ ) and ( |3l| ) can then be expressed as 
linear combinations of the one-atom trace 

Tr[p(x4-d)(x3-d)(x2-d)(xi-d)] 

where the fixed vectors x^ stand for e,e, e',e',n, or d' 
(the dipole operator of the other atom). Rather than 
calculating each term separately, we determine the one- 
atom trace for four arbitray Xq, and later substitute what 
is required by the single and double scattering terms. 



C. The single scattering vertex 

We proceed to calculate the dimensionless trace func- 
tion 

^(xa) = ]^Tr[p(x4 • d)(x3 • d)(x2 • d)(xi • d)]. (33) 

It depends linearly on the components of the Xq, , albeit 
in a complicated manner, involving the characteristics of 
the transition and the elements of the density matrix. A 
systematic way of evaluating the trace is a development 
in terms transforming under irreducible representations 
of the rotation group Q. We shall explicit the solu- 
tion in the simplest case when the atom is distributed 
with equal probability over its internal substates. Since 
the corresponding density matrix is then proportional to 
unity and therefore a scalar under rotations, the trace se- 
lects the scalar part of the averaged operator. The result 



can only be a function of the scalar products (x^ • x^j), 
of the most general form 

T(Xa) = Wl(xi • X2)(X3 • X4) + ti;2(xi • X3)(X2 • X4) 
-FW3(X1 •X4)(X2 •X3). 

(34) 

The coefficients Wi are calculated explicitly using the 
standard techniques of irreducible tensor operators (de- 
tails are given in appendix |AJ) : 



So - S2 S2-S1 S1+S2 ,„_, 

"'I = 7, I W2 = — - — , W3 = — -z — (35) 



where 



sk = 3(2Je + 1) 



1 1 K 
J J J, 



(36) 



The "6 J" -symbols Q (or Wigner couphng coefficients) 
contain all essential informations about our problem. 
They are the simplest scalar quantities that can be con- 
structed from the basic ingredients J, J^, 1 (the rank of 
the vector operator d) and K (the tensor ranks of the 
irreducible components of the scattering operator). The 
"6 J" -symbols introduce useful selection rules: 

(i) I J — Jc| < 1, the usual selection rule for a dipole 
transition; 

{ii) < K < 2: the scattering operator is the direct 
product of two vector operators and thus has irre- 
ducible components of rank K = 0,1,2. In other 
words, the change of the atomic angular momen- 
tum is limited to \m' — m\ <2 for the one photon 
scattering; 

{Hi) K < 2J: the ground state degeneracy determines 
which tensor rank comes into play. For J — 0, 
K = and thus only Rayleigh transitions m' — 
m — Q are possible; for J — 1/2, K — 0,1 and 
degenerate Raman transitions with |?7i' — ?7i,| = 1 
become possible; for J > 1, K = 0,1,2 and all 
possible transitions \m' — m\ < 2 can take place. 

A sum rule over K for the "6 J" -symbols implies that the 
w coefficients are not independent but obey 

W1+W2 + 3w3 = 1 (37) 

for arbitrary Jj_Jc. Explicit formulae for the Wi are con- 



tained in app. B 1 



We introduce a diagrammatical representation for the 
trace function (M): 




This four-point intensity vertex is the weighted sum of 
the three pairwise contractions between the vectors Xq,. 
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A factor wi comes in for the horizontal pairwise con- 
traction (xi • X2)(x3 • X4), a factor W2 for a diagonal 
pairwise contraction, and a factor ^3 for a vertical pair- 
wise contraction. It ressembles Wicks's theorem known 
from Gaussian integration |E2|, but here, the weights of 
the possible contractions arc not equal. As in quantum 
field theory, this diagrammatic representation proves es- 
pecially useful for the systematic description of higher- 
order scattering (cf. Sec. IV E). 



D. The single scattering contribution 



Using £ = xi = X4 and e' = X2 = X3 in (34), the in- 
ternal average (|2^) for the single scattering contribution 
becomes 



d^-' 



-H. 



Mj {w 



i£ 



W2\£ ■ £\ 



W3 



(39) 



The average differential cross-section for single scattering 
on an unpolarized atom is 



/das 

\ dn 



3crt, 



int 



Stt 



(Wi\£ ■£[ 



?«2 £ • £ 



W3 



(40) 



in terms of the total scattering cross-section (BJ). Us- 
ing this expression, we see that the sum rule {\i7\j simply 
represents flux conservation. All angular information is 
contained in the squared moduli of polarization contrac- 
tions. Since these expressions are even in the scatter- 
ing angle 9, the mean value (cos 9) = J d9 {da/dil) cos 9 
vanishes, justifying that the transport mean free path 
^tr = ^/(l — {cos 9)) is equal to the scattering mean free 
path i. 

To determine the bistatic coefficient now means to av- 
erage eq. (|40| ) over position. We assume a semi-infinite, 
homogenous medium of independently distributed atoms. 
The single scattering bistatic coefficient then is 



7s 



^4™ f,3/da^\ -2./ 



z>0 



\ dn U 



(41) 



The exponential takes account of the extinction of in- 
coming and scattered light with the mean free path (, 
inside the scattering medium. Since the differential cross- 
section is independent of the position, the integral is read- 
ily calculated and we find 



7s = ^ 



= -{wi\e' ■ e\'^ + W2\e' ■ e]'^ + w^) 



(42) 



The coefficients Wi{J,Jc) carry the weights of the dif- 
ferent contractions of the pola rization vectors (for de- 
tailed expressions, see app. Bl). In the case of a tran- 
sition J — 0, Jc = 1, these coefficients are simply 



{wi,W2,W3) = (1,0,0). So one recovers exactly the clas- 
sical expression 



7s 



£ •£■ 



(43) 



In the classical diagram, the upper line, read from left to 
right, signifies scattering of the wave amplitude, and the 
lower line, read from right to left, signifies scattering of 
the complex conjugate amplitude by the same scatterer 
(identified by the dashed line). The only possible connec- 
tion is horizontal, giving the factor |£' • £p that implies a 
dipole radiation pattern. For atoms, however, the coeffi- 
cients 'W2 and W3 come into play and lead to contributions 
in the lin _L lin channel (where e' •£ = £'•£ = 0) as soon 
as J > 1/2. When J > 1, there is a signal even in the 
helicity preserving backscattering channel h\\h (where 
£' • £ = 0, £' • £ = 1). Now we see why the polarizability 
( pl| ) is not sufficient to describe the scattering by a de- 
generate transition. The polarizability is essentially the 
average scattering tensor, a two-point vertex connecting 
the incoming to the outgoing polarization. If the atom 
is distributed with equal probability over its substates, 
the polarizability then is diagonal and defines a purely 
horizontal contraction proportional to Je' • £p. But the 
internal structure of an atom allows also for diagonal and 
vertical connections in the single scattering diagram: the 
classical line stretches to a two dimensional ribbon. 



E. The double scattering contributions 

Van Tiggelen et al. B^ have calculated the double- 
scattering contribution to the ladder bistatic coefficient 
in the backward direction, 
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7L2 



IQttAP 



d^rid^r2 

Zl,2>0 



-{zi+ri2 + Z2)/t 



Pl2 (44) 



'12 



for classical point scatterers in a half-space, within the 
weak scattering limit k£ ^ 1 and in the far field approx- 
imation kri2 ^ 1. Here, the exponential describes the 
attenuation of incident, intermediate and scattered light 
with mean free path £. For classical dipole point scatter- 
ers, the polarization kernel is given by Pj^2 ~ |£' • A • £p. 
For atomic scatterers under the same conditions, eq. (H) 
remains valid. As all information about the internal 
structure is connected to the polarization, only the po- 
larization kernel has to be gener alized . Keeping track of 
all factors, it follows from Sec. IV B that the polariza- 
tion kernel is given as the internal average (am over the 
square of the dimensionlcss double scattering operator: 



Pl2 = M- 



to-A-ti 



-w 



(45) 
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Using eqs. ( |29| ) to (^, it is represented by the general- 
ized ladder diagram 



Pi 



L2 



-A- 



-A- 



(46) 



This double-scattering ladder diagram is the product of 
two single scattering diagrams (Esl) connected by the po- 



larization propagator A^ 



fiiHj, one for the am- 



plitude (upper line) and one for its complex conjugate 
(lower line) . The diagram is evaluated using the following 
rules. Each scattering box yields three pairwise contrac- 
tions: horizontal with weight wi, diagonal with weight 
W2 and vertical with weight W3. Now choose WiWj for 
the two boxes and contract the vectors accordingly. For 
example, wf comes with the twofold horizontal contrac- 
tion |e • A • e'p; wiW2 and W2W1 both give |e • A • e'p. 
For the vertical connections involving factors of W3 , one 
has to use that the polarization propagator is a projec- 
tor, A • A = A, and its total contraction (arising for w§) 
is Y.^ A„; = 2. Finally 

Pl2 = (w? + wl)\e' • A • £|2 + 2wiW2\£' • A • £|2 

+ (wi + W2)w3[{e • A • e) + (e' • A • e')] + 2^3 

(47) 

For classical dipole scatterers, modeled by a transition 
J = 0, Jc = 1, one has (wi, W2, W3) = (1, 0, 0) and recov- 
ers the known result 



P, 



(cl) 



L2 



-A- 



-A- 



le' • A • £| 



(48) 



The crossed bistatic coefficient for the double- 
scattering contribution as calculated by van Tiggelen et 
al. [|43| under the same assumptions is 



7C2 



9 



IQ-kAP 



cPri(fr2 

Zl,2>0 



-(sZi+ri2+SZ2)/l 



X cos[(k + k') • ri2] Pc2 

(49) 

where s — ^(1 + 1/ cos6'). For classical dipole point scat- 
terers, the crossed and ladder polarization kernels are 



(cl) 



(cl) 



equal, Pq2 = Pl2 • This assures that in the backscat- 
tering direction (0 = 0, k' = — k), crossed and ladder in- 
tensities are equal (the strict equality for all polarizations 
is characteristic of double scattering; for higher scatter- 
ing orders, equality of crossed and ladder is only given 
for paral lel p olarizations). In the case of atomic scatter- 
ers, eq. (49) remains valid, but the polarization kernel 
has to be generalized. Casting the internal average (p^) 
in diagrammatical form, it is 



(50) 




The crossed diagram is evaluated efficiently using the 
contraction rules defined above for the ladder diagram. 
Explicitly, 

Pc2 = [wl + wl)\e' • A • ep + 2wiW3{e' ■ A • e'){e • A • e ) 
+ (wi + W3)w2[{e ■ e'){e ■ A ■ e') + {e ■ e'){e ■ A • e')] 



2wi 



(51) 



Obviously, the crossed kernel is not equal to the ladder 
kernel, Pc2 7^ Pl2- What is the relation between the 
two? In the classical theory, one habitually uses time 
reversal invariance to reduce one to the other: returning 
the lower line of the crossed diagram for classical point 
scatterers 



P, 



(cl) 



C2 



(52) 



-A-:^ 



the connecting lines are straightened out, and the crossed 
diagram becomes equal to the ladder diagram ([48| ) in 
the parallel polarization channels e' — e. This is the 
signature of reciprocity and assure a perfect interference 
contrast in the backscattering direction. But returning 
the bottom line of the generalized crossed diagram (|5C 
we find 



P 



C2 



-A- 



-A- 



(53) 



which differs from the ladder diagram (46), even if we 
put e' — e. What has happened? The ribbon that has 
replaced the classical line cannot unwind and blocks the 
diagram topologically. It blocks because the diagonal 
and the vertical contraction are not equivalent: W2 7^ it'a. 
Only in the case J = 0, Jc = 1, we have W2 = W3 — 0, 
and one recovers the correspondence to the classical point 
dipole scatterers. Eq. ( pq ) shows that W2 — w^ if and 
only if si = 0. The coefficient si stems from the anti- 
symmetric part of the scattering operator (cf. app. H). 
As t he an alysis of the double scattering amplitude in 
Sec. Ill C already showed, it is the antisymmetric part 



of the atomic scattering operator that is responsible for 
the reduction of the backscattering enhancement in the 
parallel channels. 



V. ENHANCEMENT FACTORS AND PEAK 
ANALYSIS FOR ANY ATOMIC TRANSITION 

The results of the previous section enable us to cal- 
culate analytically the intensity of polarized light scat- 
tered at first and second order by atoms which are po- 
sitioned randomly in a half-space. The spatial integrals 
in eqs. ( [44| ) and (^) are challenging because of the half- 
space geometry, but can be performed fully analytically 
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see section B3, We therefore obtain the various en- 



hancement factors at backscattering as well as the shape 
of the backscattered cone fully analytically. In this sec- 
tion, we analyze the contributions of single and double 
scattering to the backscattered intensity and determine 
the second-order enhancement factor as a function of the 
polarization channel and the atomic dipole transition. 



A. Single scattering background 

The single scattering intensityin terms of the bistatic 
coefficient 73 is given by eq. ( [42|) as a function of the 
ground state angular momentum J, the transition type 
Jo — J = 0, ±1 and the polarization vectors. Fig. Q shows 
7s for all transition types and the four standard polar- 
ization configurations. 
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0.6 

7s 0.4 

0.2 
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0.2 



12 3 4 5 

J 

FIG. 4. Single scattering intensity in terms of the bistatic 
coefficient (U3) in the backscattering direction as a function 
of the ground state angular momentum J: a) preserved and 
flipped helicity in the circular polarization channels, b) paral- 
lel and perpendicular polarization in the linear channels. Full 
symbols: parallel channels {h\\h and lin\\lin). Open sym- 
bols: perpendicular channels (hJ-h and lin 1. lin). Transition 
types: (T, V): Je = J + 1, (1,0): Jc = J, (A, A): Jo = J - 1. 

For classical point dipole scatterers, the single scatter- 
ing bistatic coefficient is 73 = | in the channels hJ-h 
and lin \\ lin (corresponding to the reflection from a mir- 
ror) and 7s = in the channels h\\h and lin ± lin [Esj. 
Fig. ^ reproduces these valu es for the transition J — 0, 
Jc — 1. As explained in Sec. IIIB| , for J > degenerate 
Raman transitions become possible and open the classi- 
cally forbidden channels: the first signal is obtained in 
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the linJ-lin channel for J = i and in the h\\h channel 
for J= 1. 

In all four polarization channels, the graphs of the two 
transition types Jo = J ± 1 (upward and downward tri- 
angles) tend towards t he sa me value as J — > cxd. Indeed, 



as shown in appendix B_l, the coefficients Wi for these 
two transition types have the same limit, corresponding 
to asymptotically equal Clebsch-Gordan coefficients. 
Two main conclusions are to be drawn from Fig. H: 

(z) A degeneracy of the atomic dipole transition leads 
to a single scattering contribution to the backscat- 
tered intensity in all four polarization channels 
(with the only exception J = ^ i^^ ^11^); this 
background signal therefore cannot be eliminated 
by polarization analysis and reduces the observable 
height (m of the coherent backscattering peak; 

(ii) The intensity in the h\\h and lin J- lin channels 
always stays below the intensity in the h J-h and 
lin II lin channels, respectively; the single scattering 
contribution thus is always minimized by choosing 
the classically forbidden channels. 



B. Double scattering interference contrast 

The contrast of second order backscattering interfer- 
ence, 



C2 



702(0) 

7L2 



(54) 



is determined by the crossed and the ladder bistatic co- 
efficients, given in eqs. (^) and ( |4g|) as integrals over 
the generalized polarization kernels ([47|) and ( pl| ) , respec- 
tively. These integrals can be evaluated analytically, and 
their expressions as functions of J and Jo in the four po- 
larization channels are contained in appendix H. Here, 
we plot the interference contrast 0-2 in Figs. ^ and |6| in 
the four standard polarization channels as a function of 
the ground state angular momentum J. 

Two features of Fig. o are particularly striking: Firstly, 
a perfect contrast C2 = 1 is obtained solely for the transi- 
tion J = 0, Jo = 1 corresponding to classical point dipole 
scatterers. The degeneracy of the atomic transition then 
degrades the contrast considerably. For instance, in the 
channel of preserved helicity {h \\ h) and for a transition 
of type Jo = J + 1 (full downward triangles in Fig. ||), 
the contrast drops to about 0.3 already at J = ^ and 
takes typical values of 0.2. Secondly, the channel h±h 
can offer a contrast up to three times higher than the 
channel h \\ h, depending on the transition type and the 
degeneracy of the ground state. 



13 



C2 



1 

0.8 
0.6 

0.4 

0.2 

1 T 

0.8 

0.6 

0.4 

0.2 






a) h\\h 



▼ ■ 

^ "^ ▼ ▼ ▼ ▼ ▼ T ♦ 

A^ ,^ ,. ^ ^ ^ ^ '^ '^ 



V V 



A^ Q-Q 



h) h±h 
^ -V -V V I 

A- A A -^-i 



2 3 

J 



FIG. 5. Contrast of double backscattering interference (p4 
as function of the ground state angular momentum J for cir- 
cular polarizations: a) conserved helicity, b) flipped helicity. 
Transition types: (T, V): Jc = J + 1, (■,□): Jo = J, (A, A): 

Je = J - 1. 
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FIG. 6. Contrast of double backscattering interference (m) 
as function of the ground state angular momentum J for lin- 
ear polarizations: a) parallel polarizations, b) perpendicular 
polarizations. Transition types: (▼, V): Jc ~ J + 1, (■.□): 
Jo = J, (A, A): Jc = J- 1. 

Figs, ^d) and g|b) show that in the crossed channels 



(/i _L /i and lin _L lin) , the contrast is always maximized 
for a transition type Jo — J -f 1. But in the parallel 
channels {h\\h and lin \\ lin) and for larger values of J, 
the contrast is optimized for Jo — J. In the limit Jo — 
J —^ oo, the contrast C2 even approaches 1. A contrast of 
1 indicates that the antisymmetric part of the scattering 
tensor vanishes. Indeed, the Clebsch-Gordan coefficients 
display a symmetry that suppresses the antisymmetric 
part of the scattering tensor as Jo = J — > oo. 



C. Backscattering enhancement factor 

Fig. H shows that the smallest single scattering signal 
is obtained in the h\\h channel for Jo = J -I- 1. This 
configuration could also be expected to render the best 
enhancement factor. However, Fig. ^ shows that the in- 
terference contrast in this configuration is particularly 
low. As will indeed be seen in this section, the choice of 
Jo = J + 1 and h\\h does not guarantee an optimized 
backscattering enhancement. Depending on the degener- 
acy, the crossed channel or the transition type Jc — J can 
offer a better interference contrast and lead to a higher 
enhancement factor. 

An enhancement factor up to second order, 



a = 1 



7C2 (0) 

7s + 7L2 ' 



(55) 



combines the single and double scattering contributions. 
It has to be pointed out, however, that its exact value 
may not be compared directly to experimental results. 
Indeed, either the scattering medium has the semi- 
infinite geometry of a half-space, but then third and 
higher scattering orders cannot be neglected. Or it has 
the finite geometry of a laser-cooled atomic cloud which 
truncates higher scattering orders, but then the relative 
weight between single and double scattering is modi- 
fied. Numerical simulations can determine the role of re- 
stricted geometry and are currently under study. Prelim- 
inary results show that the ratio of the double-scattering 
crossed intensity to the ladder intensity is almost inde- 
pendent of the shape of the medium. In other words, it is 
the internal structure which is essential for the low con- 
trast of the interferences, not the spatial arrangement of 
the various atoms. Thus, the present analytical calcula- 
tion permits to follow how the effects of single scattering 
background and reduction of interference contrast com- 
bine to result in small enhancement factors. 

Figs. u\ and B show the enhancement factor (55) as a 
function of the ground state angular momentum for the 
four standard polarization channels. 

The difference between Figs, g and is given by the 
single scattering contribution, shown in Fig. Ha). In the 
channel of conserved helicity {h\\h), the lowest single 
scattering intensity is observed for the Jo = J -I- 1 tran- 
sition type, so that the enhancement factor, following 
closely the interference contrast, drops from its classical 
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value 2 to about 1.2. The already poor contrast for the 
Je = J — 1 transition type is further reduced by single 
scattering. The increasingly good contrast for Jc = J 
at higher values of J is counterbalanced by an impor- 
tant single scattering contribution, so that the effective 
enhancement stays below 1.4. We thus find that the clas- 
sical enhancement factor of 2 in the helicity preserving 
channel is irrevocably lost for atomic scatterers as soon 
as J > 0. 

Although the contrast of interference tends to be 
higher in the hJ-h channel than in the h\\h channel 
(Fig. ^, the single scattering contribution (Fig.^b)) also 
is more important, resulting in low enhancement factors 
below 1.31. 




FIG. 7. The second-order backscattering enhancement fac- 
tor (B5f) as function of the ground state angular momentum 
J for circular polarizations: a) conserved helicity, b) flipped 
helicity. Transition types: (T, V): Jc = J + 1, (■,□): Jc = J, 
(A, A): Jc = J-1. 

The enhancement factors in the linear channels, dis- 
played in Fig N, show the same characteristics. With the 
only exception of the J = 0, Je = 1 transition in the 
linl. lin channel, we find that all possible atomic transi- 
tions yield enhancement factors below 1.35. 

The interplay between single scattering background 
and interference contrast makes it difficult to predict in 
which configuration the optimal enhanced backscattering 
can be measured. Intuition formed with classical scatter- 
ers would recommend a transition of type Je = J -|- 1 and 
the h\\h channel. But for a high enough degeneracy of 
the atomic transition, classical intuition turns out to be 
a bad counsellor. For J = 3, Je = 4, the calculated 
effective enhancement factor is higher in the h J-h chan- 
nel (a = 1.21) than in the h\\ h channel (a = 1.17). This 



had first been observed experimentally jlj] and remained 
puzzling until taking into account the atomic internal 
structure [031. 
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FIG. 8. The second-order backscattering enhancement fac- 
tor ( pq ) as function of the ground state angular momentum J 
for linear polarizations: a) parallel polarizations, b) perpen- 
dicular polarizations. Transition types: (T,V): Je = J-|-l, 
(■,n): Jc = J, (a, a): Jc = J - 1. 

Fig. M indicates that an optimized enhancement is ex- 
pected for a transition Jc = J in the h\\h channel. How- 
ever, a direct experimental verification seems delicate be- 
cause a transition of type Jc = J is not closed in general 
(the emission of a photon from the excited level Jc to 
a final level Jc = J — 1 is allowed). These events cut 
off elastic scattering paths and yield a high background 
intensity, unfavourable for experimental detection. An 
interesting exception to this rule is the closed transition 
J = Je = 21 which has the additional advantage that no 
single scattering background pollutes the h\\h channel. 



D. Enhanced backscattering peaks for J = 3, Jo = 4 



The scattered intensity enhancement 



a(/i) = 1 



7C2(m) 

7S + 7L2 



(56) 



as a function of the reduced scattering angle fi = ki6 for 
any atomic transition and any polarization is given ana- 
lytically in terms of the bistatic coefficients (see app. ^ 
for details). Fig. H displays the backscattering peak a(/i) 
for the case that has been experimentally studied: the op- 
tical transition between two hyperfine levels (F = 3 and 
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Fe = 4) of laser cooled Rubidium atoms jlj] . Hyperfine 
levels are characterized by a total angular momentum F 
including the coupling with the nuclear spin; our anal- 
ysis applies to any total angular momentum which we 
continue to note here by J. The highest peak arises in 
the channel of flipped helicity {hJ-h), the linear peaks 
are almost equivalent, and the smallest peak is given for 
preserved helicity {h\\h). The calculated enhancement 
factor in all four channels is of the order of 1.2. The ex- 
perimentally measured enhancement factors are smaller 
than the present values because the atomic cloud has nei- 
ther a uniform density nor the geometry of a half-space. 
Nevertheless, we stress that the calculated peaks repro- 
duce semi-quantitatively the experimental ones as shown 
in||]. 

The shape of the CBS cone for atoms is similar to the 
one for point dipole scatterers: the angular width is of the 
order of l/k£ and can vary by a factor 2 depending on the 
polarization channel. In the helicity channels, the CBS 
cone is isotropic. In the linear channels, its presents an 
anisotropy which is characteristic of polarization mem- 
ory in low-order scattering [Q, see also App. B 5| . In 
the tin \\ tin channel, it extends further in the direction 
of the polarization than perpendicularly, reflecting the 
anisotropy of the dipole scattering cross-section in the 
Fourier plane. 
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FIG. 9. The intensity enhancement (|56| ) as a function of 
the reduced scattering angle /i = k£0 for J = 3, Jo = 4 in 
the four polarization channels. In the linear channels, the 
intensity is scanned in the direction parallel to the incident 
polarization. In the helicity channels, the intensity is inde- 
pendent of the scan direction. The dashed curve in the h\\h 
channel shows CBS contribution from Rayleigh transitions 
{m' = m) only. In this channel, the dominant contribution to 
the CBS peak comes from Raman t ransitio ns between differ- 
ent substates m' 7^ m, see also Sec. [II C 3. 



tering paths dominate for larger scattering angles. An- 
alytical results for double scattering thus provide infor- 
mation about the wings of the backscattering peak that 
are in principle measurable experimentally. 



App. 



B5 



contains the analytical expressions for the wings of the 
backscattering peak in the four usual polarization chan- 
nels. It can be seen that the intensity decreases as 
{k£9)~^ in all four channels, with coefficients depend- 
ing on J, Jc. Furthermore, the anisotropy in the linear 
polarization channels, i.e., the dependence of the scat- 
tered intensity on the angle between incident polarization 
and the direction of the scan, decreases as the degener- 
acy of the atomic transition increases. This reduction 
of the anisotropy is consistent with the intuitive picture 
that degenerate atomic transitions depolarize the inci- 
dent light more efficiently than dipole point scatterers 
and that the memory of direction as seen in the backscat- 
tering anisotropy is lost more rapidly. 



VI. CONCLUSIONS AND OUTLOOK 

The internal structure of an atom determines its light 
scattering properties. We have shown that a degeneracy 
of the atomic dipole transition reduces the observable 
backscattering enhancement factor in two ways: single 
scattering is present in all polarization channels and the 
non-scalar part of the scattering tensor reduces the con- 
trast of CBS interference. A complete analytical solution 
for the case of unpolarized atoms has been presented to- 
gether with a generalization of the classical ladder and 
crossed diagrams to the case of atoms. An immediate 
extension, under current study, is the application of the 
present method to higher orders of scattering. Finally, 
going beyond the weak localization regime, further re- 
search is needed in order to decide whether the the inter- 
nal structure is not a substantial difficulty in the quest 
for the strong localization of light in cold atomic gases. 
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A small scattering angle is associated to endpoints of 
scattering paths lying far apart. Conversely, short scat- 
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APPENDIX A: TRACE EVALUATION USING IRREDUCIBLE TENSOR OPERATORS 

In the following, we employ the standard theory of irreducible tensor operators as exposed in the textbooks by 
Edmonds [^ and Blum [^. We have to calculate 

T(xo) = ]jJ;Tr[pj(x4 • d)(x3 • d)(x2 • d)(xi • d)] (Al) 

where d = d'^^ is the reduced dipole operator d = D/D, an irreducible tensor operator of rank 1 acting upon 
the eigenstates \Jm) of the angular momentum operators J^ and Jz. Its reduced matrix element is by definition 
(Je||d||J) = \/2Jc + 1- We introduce the ratio of multiplicities Mj — (2Jc + 1)/3(2J + 1) for convenience. The 
Xq = Xq are fixed free vectors commuting with J^, Jz,d. Therefore, the Xq are irreducible tensors of rank 1, but 



not operators, and the trace (Al) acts only on d. 



Let O be an operator decomposed into its irreducible components, 

= ^aL,0(^). (A2) 

L,q 

Its average (O) = TrpO in a system described by a density matrix p can again be decomposed, 

TrpO = 5]-^^5]p^(J,J')(J'||0(^)||J). (A3) 

L,q ^ J,,}' 

All angular information has been concentrated into the coefficients aLq and the components 

PqiJ,J')^ ^(-)^-'^'-™(JJ'-mm'|L9)(Jm|p|J'm') (A4) 

of the so-called statistical tensor operator. 

If the system is distributed with equal probability over all substates \Jm) for a given J, the density matrix with 
elements {Jm\p\J'm') — (2J + ^)~^Sjj>dmm' is purely scalar, and its only non-zero irreducible component is p^ = 



(2J + l)-i/2^ The trace (^ then reduces to 



TrpO = -^^^(J||0(°)||J). (A5) 

All we have to do now is to decompose the operator O = (x4 • d) (xa • d) (x2 • d) (xi • d) into its irreducible components, 
determine the coefficient aoo and the reduced matrix element (JHO'"^!! J). We begin with decomposing the operator 
of second order 021 = (x2 • d)(xi • d). The scalar products can be expressed in any basis, in the Cartesian basis as 
well as in irreducible components. 



K 

021 



= E E (-)^-"[-2Xi]L^jMd]L^). (A6) 



K=Om=-K 



Here, [A'-'^^B'-'' '] denotes the direct product of two irreducible tensors. The irreducible components of the product 
are composed from the irreducible components of the factors, 

j^(fe)^(fe')](/0 = Y, {kk'rs\Km) A^.'^^ Bp , (A7) 

r,s 

using the Clebsch-Gordan coefficients {kk'rs\Km). Application of the inverse formula 

A[^'B^S^ = ^ {KK'mm'\Lq) [A^^^ B^'^X'^ (A8) 

L,q 



to the product O = 043O21 leads to the decomposition 



[L) 



q 



K,K',L,q ^"^ 
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which is a hnear combination (sum over K,K') of totally decomposed operators (sum over L.q). Under the trace 



according to eq. (A5), only L = q = survives so that we are left with a sum of three terms K = K' = 0,1,2. 



The reduced matrix element {J\\0'^°\K)\\J) = (J||[[dd]('^)[d(i](^)](°)|| J) can be calculated using the general formula 

(j'||[^(fe)5(fe')]('c")||j) ^ {-f'+-'+J'{2k" + l)^^^Y.{'j J' J" I {J'W^''^\\J"){J"\\b'^'''^\\J) (AlO) 



for the reduced matrix element of the direct product of two irreducible tensor operators acting on the same system [ pO[ . 
Two iterated applications of this formula yield 

(J||0'«)(A-)||J).(2J. + l)V)''-(f^)"'{,l i 'If^ (All) 

The last thing to do now is to evaluate 

i~f{2K + iy/\oo{K) = [a;iX2](^) • [x^x^]^''^ (A12) 

This scalar product of two irreducible tensors of rank 2 can again be written in any basis, in irreducible components 
as well as in Cartesian components, 



[XiX2 



(^).[x3X4](^)=5:[xi.2]!f[a.3:.4]lf. (A13) 



iW 



The cartesian components [xQa;/3]^ are given by the usual decomposition of matrices: for K = 0, the scalar part or 
trace 

[XaXf3]fj'' = -(X„ • X0)%, (A14) 

for K — 1 the antisymmetric part 

[XaXpl^j ^ '^{XaiXf3j XajXfji) (Aioj 

and for K = 2 the traceless symmetric part 






[XaXp]}/ = -{XaiXfjj + X^jXpi) - -(Xq • X^)%. (A16) 



Putting everything together, we summarize 

2 



r(x,) = 3(2Je + i)^|) ) ;|| Tk(xO 



K 

To = -(Xi-X2)(x3-X4) (A17) 

Tl = 2 t'-^l ' ^4)(X2 • X3) - (Xi • X3)(X2 • X4)] 

T2 = 2 [(^1 • X4)(X2 • 3:3) + (Xl • X3)(X2 • X4)] - Tq. 

This form shows nicely that the scalar, antisymmetric and traceless symmetric parts of the single scattering operator 
combine with their counterparts in the direct product and contribute to the scalar trace. Regrouping of the different 



contractions leads to the vertex form (p4) presented in Sec. IV B 
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APPENDIX B: ANALYTICAL EXPRESSIONS OF 
THE DOUBLE SCATTERING CONTRIBUTIONS 

1. Values of transition-dependent coefficients 

For Jo = J + 1, the transition-dependent coefReients 
(pq) are explicitly 



1 



io(j + i)(2j + i; 



6J2 + 17J+10, i = l 

X { -4J(J + 2), i = 2 

J(6J + 7), ? = 3 

(Bl) 



For Jc = J, 



1 



2J2 + 2J+1, 



i^ 1 



^ ^ '2J2 + 2J+1, i = 3 



For Jp = J- 1, 



1 



10J(2J+1) 



(6J+1)(J-1), z = l 
-4(J+1)(J-1), 1 = 2 (B3) 
(J + 1)(6J-1), i = 3 



As pointed out in Sec. IV E , the antisymmetric part of the 
scattering tensor plays no role when W2 — W3. The only 
finite values of J, Je for which this condition is fulfilled 
are J = 0,Jc — I, the case of the classical dipole point 
scatterer, where 'W2 = ws = 0. The coefficients take non- 
trivial values in the limit J -^ 00: 

1 [(3,-2,3) ,Je = J±f ,„,, 
^ ' ' 10 [(2,2,2) ,Jc = J ^ ' 

and we see that a non-trivial realization W2 — W3 — ■^ oi 
a vanishing antisymmetric part of the scattering tensor 
is given asymptotically in the case Jc = J — > 00. 



2. Ladder contribution 

The six-dimensional integral M4) with the generalized 
ladder polarization kernel (^) can be exactly calculated. 
The first (trivial) step is to use the translational invari- 
ance perpendicularly to the incoming direction, and re- 
duce it to an integral over the three components of the 
interparticle vector ri2 and over (zi -I- Z2)/2. In a sec- 
ond step, we use spherical coordinates (ri2, t9, (p) for ri2, 
where •& is the angle between the 2;-direction and ri2, and 
if the azimuthal angle (in the circularly polarized case, 
the ladder kernel is independent on ip). The integrals over 
ri2 and (zi -I- 22)/'^ are then easily performed, leading to 
the double scattering ladder contribution: 



7L2 



9 
32^ 



sinz? Pl2(^,<^) dMcp 
1 + |cosi?| 



(B5) 



expressed as an integral over the direction n — 
(sin i9 cos </?, sin sin (/3, cos ^) of the interparticle vector. 
The kernel Pl2, given by eq. (p7|), involves only simple 
trigonometric functions of i^ and ip, which makes the cal- 
culation of the integral easy. The result depends of course 
on the incoming and outgoing polarizations e and e' . 
We finally obtain: 

9 

7l2 = ^ ('1(^1 + W2)'^ + I2W1W2 + ^3(^1 + W2)'Wj, + hwl) 

(B6) 



Here, the terms 



have been completed to (wi 



1102)'^, simplifying all following expressions. The coeffi- 
cients Wi{J, Je) carry the dependence on the atomic tran- 
sition, and the coefficients k are given as functions of the 
polarization channels: 



(B7) 





h\\h 


h±h 


lin lin 


lin _L lin 


^1 
12 


5 

48 

21n2-l 


ln2- i^ 
-(21n2-l) 


1^2 -i 



5 

96 




and 



Za = 2 In 2 - 



1 



/4 = 2 In 2 



(B8) 



in all four channels. The coefficient li for the four chan- 
nels had been derived in ref . |Q , the others describe the 
generalization to the case of degenerate atomic transi- 
tions. 



3. Crossed contribution 

The calculation of the crossed contribution given by 
eqs. (E^) and (pl[) follows the same lines. There is how- 
ever a complication due to the cos[(k + k') • Y12] term. 
We choose the spherical coordinates (7-12, ??, </?) such that 
the X axis is along (k + k'), that is in the direction of ob- 
servation. The integral over the transverse components 
of ri2 and over (zi -I- 2:2)72 and ri2 yields the following 
result: 



7C2(/i) = — - 
OZTT 



where 



sini9 Pc2{"d,p) dMcp 



1 + I C0S1?| -|- /U^(l — I C0S??|)C0S^ if 



M 



eke 



(B9) 



(BIO) 



is the reduced scattering angle. 

The kernel Pc2 is a combination of simple trigonomet- 
ric functions of 1!) and ip. This makes it possible to cal- 
culate easily the integral over p, leading to: 



7C2(Ai) 



1 
dx 



Cy(X', J, Jo) 



^(H-a;)2+^2(i_2.2) 



(Bll) 
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where the crossed kernel C{x;J,Jc) depends on the 
atomic transition J —> J^ via the coefficients w^, 



C(x; J, Jo) = {Wi + W3)^Ci{x) + WiW3C2{x) 
+ {Wi + W3)w2C3{x) + wlc4{x) 



(B12) 



and the functions Ci{x) depend on the polarization chan- 
nel: 





h\\h 


hJ-h 


lin lin 


lin -L lin 


ci(x) 


ki^-xr 


\{l + xr 


1(1+^2)2+^11 


A± 


C2ix) 


2x^ 








2x^ 


csix) 


l + x"^ 





1 + x^ + B\\ 





ca{x) 


2 





2 













(B13) 



In the h J-h channel, the only non-zero coefficient is 
ci (x) . This means that apart from a multiplicative factor 
(wi + ws)^, the backscattering peak for any atomic tran- 
sition has exactly the same shape as the classical peak. 
This is due to the fact that only Rayleigh transitions con- 
tribute to the CBS peak in the hJ-h channel and that the 
radiation diagram of such transitions is - averaged over 
the magnetic quantum number - identical to the one of 
classical point dipole scatterers. In all other channels, the 
form of the backscattering peak itself is changed, be it 
only in minor ways. In the linear channels, a supplemen- 
tary complication arises because the intensity depends on 
the angle (f> between the incident polarization vector and 
the direction of the observation. </> = corresponds to a 
scan parallel to the incident polarization vector, cj) — ^ 
to a scan perpendicular to the incident polarization vec- 
tor (in ref . ^] , the opposite convention is chosen) . This 
anisotropy of the backscattering enhancement, observed 
already for classical point scatterers, is contained in the 
expressions 



(1 



^2^2 



-(l+X^COS' 



All = 



B« = {1 - x'^)X cos2(j) 



1 



-X cos 20 



(B14) 



A, 



(1 



^2\2 



-(l-X^COS' 



with 



y_, ^ ^{l + x)^+^^l-x^)-l-x 
^-^-^ ^3^p • (B15) 



Finally, the integral (Bll) can be calculated analyt- 
ically. The expressions are rather complicated and we 
give them for completeness. We obtain: 



7C2(m) = (^"1 + W3)^7i(m) + WlW372(Af) 

+ {wi + W3)w2j3{n) + wl-fiiy) 



(B16) 



where the non-zero functions 7i(/i) are given by the fol- 
lowing expressions: 



In the h\\h channel: 
3 



71 (m) 

72 (Ai) 

73 (m) 



256(1-^^2 



[32 - 176^*^ - 84/ + 18fi^ 



+ (-22 -I- 144/ - 17/)v/l + / 
+ 3/(48 - 16/ + 3/)F(^)] 
9 



8(1 -/)2 
+ (2 + /)F(/] 
9 



[-4- 2/ + 3^1 + / 



16(1 -/)2 
+ (4 - 4/ + 3/)i^(M)] 



[-4-2/ + 3Vr+/ 



74 (a^) = -^F{^i) 



• In the hJ-h channel 
3 



(BIT) 



71 (/ 



J [ - 2(80 - 56/ + 42/ + 39/) 



71 (/ 



256(1-/) 
+ (122 - 144/ + 127/)v/l + / 
+ 3(32 - 64/ + 96/ - 48/ + 19/)F(/] 

(B18) 

In the lin || lin channel: 
3 



[ - 288 + 48/ - 252/ - 138/ 



73 (/ 



512(1-/ 

-1- (222 - 144/ + 237/) v/TT/ 

+ (192 - 384/ -I- 720/ - 336/ + l23fi^)F{ii) 

+ Ai (^) cos 2(j) + A2 ifi) cos 40] 

:^^^^[-4-2/ + 3v/TT/ 

+ (4 - 4/ + 3/)F(/i) + Biii) cos 2</)l 



74(/ = -^Fiji) 



• In the lin _L lin channel: 
3 



(B19) 



71 (m) 



[32 - 176^*^ - 84/ + 18/ 



72 (a*) = 



512(1 -/)4 
+ (-22 + 144/ - 17/) VTT/ 
+ 3/(48 - 16/ + 3/)F(^) 
~A2in) cos 4:(j)] 
9 



(l-/)2 

(2 + /)F(/] 



[-4-2/ + 3v/l + / 



(B20) 



All other 7^ are zero as evident from (B13). In the l inear 
channels, the anisotropic contributions from (B14) are 
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weighted by 

Aiifi) = [ - 56(-2 + 8/1^ + 4/i'' + 5//) 

+ 28(-4 + IS^i' - 14/ + 15Ai^)v/l + /i2 
+ 12//(16 + 8/i' + 6// + 5/i'^)i^(Ai)] /m' 
A2(/i) = [48 - 152^2 + 128/ + 48//^ - 212/ - 70^^° 

+ (-48 + 176/ - 222/ + 88/ + 111/) VTT/ 
+ 3/(8 + 24/ + 3n^)F{^i)] // 
B(^) = [2 - 4/ - 4/ + (-2 + 5/) VTT/ 



-/(2 + /)F(/]// 



(B21) 



In all these expressions, the auxiliary function F{jj,) is 
given by: 



F(n) — 2argcosh ( i—r ) — argcosh ( ^r 



(B22) 



Under this form, F{p) is not a manifestly real function 
of /i. It can be rewritten as: 



7T=;^^'' [ V2/ ^^^ 



|Ai|<l 



IVm^^ 



2 / Vl + / - 1 ^^^ \ , , 

arcsinl-^^ ^- V/ - 1 , Ia^I > 1 



V2^2 



(B23) 



In table (B13), the function C4(x) is just a constant, with- 
out any angular dependence on (/) o r a; — cos??; The 

corresponding contribution 74(/i) in (p317) and (B19) is 



essentially F{^). We see therefore that F{^) is - within a 
factor 9/4 - the crossed double-scattering bistatic coeffi- 
cient for a scalar wave scattered by a semi-infinite homo- 
geneous medium of point scatterers. It is a bell-shaped 
function around fi = with width of the order of unity. 
With the help of the previous expressions, the scat- 
tered intensity can be plotted, for all atomic transitions, 
all polarization channels and all directions of observation. 



4. Crossed contribution for exact backscattering 

In exactly the backscattering direction /i = 0, the 
above expressions simplify considerably, yielding the 
crossed bistatic coefficient in the backscattering direction 



9, 



7C2(0) = ^(ci(u)i + W3)' + C2W1W3 
o 

-I- C3(wi + Wz)w2 + C4W2) 



(B24) 



where the numerical coefficients Ci are given as functions 
of the polarization channels: 





h\\h 


h±h 


lin lin 


lin -L lin 


Cl 
C2 
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48 

21n2-l 


ln2-i 



ln2-i 



5 

96 

21n2-l 


C3 


21n2-i 





21n2-i 





C4 


21n2 





21n2 






Just as for the ladder contribution, the coefficient ci had 
been derived in ref. Q], the others describe the gener- 
alization to the case of degenerate atomic transitions. 
In the parallel channels, all crossed coefficients c,; are 



equal to the corresponding ladder coefficients U (B7-B8). 
This is the signature of reciprocity since the ladder and 
crossed contributions are then equal for W2 — W3. In the 
perpendicular channels, no such correspondence can be 
observed. 



5. Wings of the crossed contribution 

For a large reduced scattering angle fj, — k£9 ^ 1, the 
previous expressions can be expanded in powers of /i~^, 
giving the wings of the enhanced backscattering peak. 
This asymptotic expression describes the wings of the 
backscattering peak even if higher orders of scattering 
contribute to the intensity at smaller angles. The crossed 
bistatic coefficient in the wings becomes 



97r 
7c (m) = s-iaiiwi + W3)' + a2WiW3 
o/i 

-I- a3{wi + W3)w2 + aiwl) + 0{^i~'^) 
where the wing coefficients are 



(B25) 
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i(l+C0s2,/)) 





04 
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The wing coefficients in the linear channels depend on 
the angle (j) between the incident polarization and the di- 
rection of the intensity scan, carrying the anisotropy of 
the linear backscattering peaks. In the lin \\ lin channel, 
the intensity is higher in the direction of the polarization 
{(j) — 0) than perpendicular to it, yielding a cigar-shaped 
intensity pattern in the plane of observation (in eq. (A. 5) 
of ref. 1^^, a term — |sin^20 is missing, otherwise all 
coefficients ai coincide). In the lin J- lin channel, the in- 
tensity is smaller in the directions of incident {4> = 0) 
and scattered {(j) = tt/2) polarization than along the di- 
agonals, yielding a cloverleaf pattern in the observation 
plane. As pointed out in ref. p3| , in the lin J- lin chan- 
nel the classical coefficient ai oc sin^ 2<j) vanishes if the 
intensity is scanned in the direction parallel or perpen- 



dicular to the incident polarization 



0, f). That 



means that the peak decreases in these directions as /i 
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instead of /i ^. But for atoms, a second constant coef- 
ficient a2 — h comes into play that maintains a (mod- 
ulated) decrease as /i^^ in all directions, thus reducing 
the anisotropy In the lin \\ lin channel, the above expres- 
sions permit to verify that the classical anisotropy ratio 
7c (0 = 0)/7c((/> = f ) = I decreases as the atomic de- 
generacy increases, converging to H for transitions of 



type Jo 
as J -^ oo. 



J ± 1 and to 22 for transitions of type Jc = J 
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